Orbital-Peierls State in NaTiSi 2 O e 
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Does the quasi one-dimensional titanium pyroxene NaTiSi2 06 exhibit the novel orbital-Peieris 
state? We calculate its groundstate properties by three methods: Monte Carlo simulations, a spin- 
orbital decoupling scheme and a mapping onto a classical model. The results show univocally that 
for the spin and orbital ordering to occur at the same temperature -an experimental observation- 
the crystal field needs to be small and the orbitals are active. We also find that quantum fluctuations 
in the spin-orbital sector drive the transition, explaining why canonical bandstructure methods fail 
to find it. The conclusion that NaTiSi2 0e shows an orbital- Peierls transition is therefore inevitable. 
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Introduction. Transition-metal compounds display a 
dazzling collection of physical properties, spanning the 
range from colossal magneto-resistance in manganites, 
the spin-Peierls effect in quasi one-dimensional systems 
to high T c superconductivity in cuprates 0, • In these 
materials the interactions between electrons can be so 
strong that they destroy the Fermi surface. A Mott insu- 
lator is the canonical example. In this situation electrons 
become localized and they regain their local spin, charge 
and orbital degrees of freedom. The orbital degrees of 
freedom play a special role here since they can couple 
both to the lattice, via the Jahn- Teller effect, and to 
the spins, via superexchange interactions [3]. The inter- 
play between all these degrees of freedom is the apparent 
root of the observed wealth in symmetry broken ground- 
states in e.g. doped manganites and even the driving 
force behind multiply ordered systems, such as e.g. mul- 
tiferroics |j,|5|- 

Recently a novel type of groundstate, the orbital- 
Peieris state, was proposed to exist in several compounds. 
This state is expected to appear in one-dimensional sys- 
tems where orbital degrees of freedom lead to the dimcr- 
ization of the atoms that make up the chain, in anal- 
ogy to the spin-Peierls state, where the dimerization is 
driven by the electron spins. As the orbital degrees of 
freedom couple to both the lattice and spin, the orbital- 
Peieris state is expected to reveal itself experimentally 
through orbital, spin and lattice anomalies. The search 
for this state, specifically in vanadates and titanates, mo- 
tivated a number of experimental and theoretical stud- 
ies H 0, H ©, H El 13 13 13 111 13 E3 and the 
findings became the subject of intense debate. The 
first candidate with possibly an orbital-Peieris ground- 
state is the perovskite vanadate YVO3. However, the 
evidence in favor H Q is opposed by extended band- 
structure calculations that indicate that the compound 
is electronically three-dimensional ||. In this letter we 
focus on the other possible candidate 0|j the pyrox- 
ene compound NaTiSi2C>6. Studies of its microscopic 



spin-orbital Hamiltonian have shown that indeed this 
novel state can form, but ab initio LDA bandstruc- 
ture and ionic cr yst al field calculations contradict this 
claim I3I3H13QI3I3G3- Here we present a 
unified picture and show that indeed it exhibits orbital- 
Peieris ordering. 

The titanium pyroxene. The distinguishing feature of 
the pyroxene compound NaTiSi2C>6 is the presence of 
quasi one-dimensional arrays of edge sharing TiOg oc- 
tahedra. Besides its spin (s=l/2), each Ti 3+ ion has 
an additional orbital degree of freedom. The system 
shows a dramatic drop in the magnetic susceptibility be- 
low Tqo = 210 K [lfj, while at the same temperature 
a structural change is observed. Raman scattering re- 
veals that this transition is accompanied by a phonon fre- 
quency shift and broadening [111 Il4| . It was immediately 
suggested that this apparent connection between lattice 
and spin degrees of freedom is driven by the ordering of 
titanium ti g orbitals. Starting from a tight binding Hub- 
bard model a microscopic spin-orbital model was derived 
in support of this picture and it was shown explicitly 
that the spin-orbital model can have an orbital-Peieris 
groundstate at low temperatures 0, 0] . 

However, also two pieces of evidence disputing this 
claim have appeared. Ionic crystal field calculations sug- 
gest that the ti g orbital splitting is substantial in this 
compound Since the microscopic models are based 
on the assumption that the orbitals are (nearly) degen- 
erate, these results seem to strike a serious blow to the 
proposed orbital-Peieris picture. Besides this, density 
functional calculations by Popovic et at indicate that 
a Haldane spin-one chain forms at low temperature, pre- 
empting the orbital-Peierls transition Recently this 
result was challenged by LDA+U bandstructure calcu- 
lations, where electron correlation effects are partially 
incorporated on a mean field level [l7j . 

In this Letter we will show that we can account for 
all of these results, by an analysis that is based on the 
microscopic model that one of us originally proposed in 
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Ref. ^lj. Both numerical Monte Carlo simulations and 
an analytical spin-orbital decoupling for the microscopic 
model show that even if the crystal field splitting is com- 
parable to the size of the inter-site exchange coupling, 
the orbital-Peierls transition still occurs. Thus for the 
spin and orbital ordering to occur at the same tempera- 
ture -as it does experimentally- the crystal field should 
be smaller than or at most comparable to the inter-site 
exchange energy. This places a strict upper bound on the 
crystal-field splitting. Furthermore we find that quantum 
fluctuations in the spin and orbital sector actually drive 
the orbital-Peierls transition. Precisely these correlations 
are not taken into account in canonical band-structure 
calculations. Neglecting them in our microscopic Hamil- 
tonian actually reduces it to a solvable classical model. 
An exact transfer matrix computation shows that with- 
out the quantum fluctuations indeed alternating ferro- 
magnetic correlations between neighborin g sp ins appear, 
as in bandstructure calculations of Ref. 13] . This ex- 
plains the failure of this method to find the orbital-Peierls 
groundstate. Judging from these results the conclusion 
that NaTiSi2C>6 exhibits an orbital-Peierls transition is 
inevitable. 

Microscopic Model. Each Ti 3+ ion in the edge-sharing 
TiOe octahedra has a single electron in one of the three 
tig orbitals denoted by \xy), \yz) and \zx). The \zx) 
orbital turns out to be inert, which renders it higher in 
energy and it is thus neglected. The Coulomb interac- 
tion between electrons on the same titanium atom is so 
large that the exchange interactions can be determined 
by a second order perturbation expansion in the elec- 
tron hopping parameter. The leading order spin-orbital 
Hamiltonian is |TT| : 
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where we use spin operators S and orbital operators T 
(where T z = ±i and the plus (minus) sign correspond 
to an occupied \xy) (\yz)) orbital). Nearest neighboring 
lattice sites are denoted by i and j and J denotes the 
exchange integral. Here we neglected the small contribu- 
tions due to the Hund's rule exchange, but we checked 
for each case that our final results do not depend on this 
approximation. We also include two orbital-only contri- 
butions to the Hamiltonian, the crystal field splitting, 
lifting the degeneracy of the two ti g orbitals, parameter- 
ized by Jcf, and a term due to the superexchange [l2j |: 
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Methods. We study the model given by equations i ffl l 
and J2J) by three different methods. First we performed a 
Monte Carlo simulation of the Hamiltonian for a system 
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FIG. 1: Left: orbital and spin transition temperatures ver- 
sus crystal field splitting determined from the Monte Carlo 
calculations. In the region left of the dotted line the spin 
and ordering temperature coincide. Right: transition tem- 
peratures versus crystal field splitting from the spin-orbital 
decoupling calculations. 



of 100 sites. In this simulation we added an inter-chain 
coupling in a mean field fashion, in order to account for 
the three dimensionality of the system: 



H = H 
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(T t z ot ). (3) 



Here Jinter is the effective mean field inter-chain cou- 
pling, z is the number of neighbor sites and J' the bare 
inter-chain coupling strength. We take J' = J/10 in the 
present calculations. It is worthwhile to note that one 
can use a classical Monte Caro simulation for this Hamil- 
tonian, rather than a quantum Monte Carlo algorithm. 
This is justified by a particular symmetry in the spin- 
orbital model: the staircase structure of the titanium 
chains combined with the orientation of the t? orbitals 
does not allow any spin states other than single spins or 
dimcrs. 

Within the Monte Carlo simulation we determine the 
magnetic (spin) and orbital susceptibilities, \ S an d X° ■> 
as well as the nearest neighbor correlators for spins and 
orbitals: C s and C° respectively. From the tempera- 
ture dependence of the spin/orbital susceptibility we can 
then also deduce the transition temperatures below which 
spin and orbital ordering are found {T s and T°). These 
quantities we studied as a function of both J and Jcf • 

In the second method we employ a decoupling of the 
spin and orbital degrees of freedom of the Hamilto- 
nian QJ. We decouple Hst by introducing static orbital 
and spin fields: 
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FIG. 3: A sketch of the zero temperature states. The lobes 
of the orbitals that lie in the plane of the titanium chain are 
depicted. The dark lobes are occupied orbitals. Top: orbital- 
Peierls state with spin valence bonds. Bottom: alternating 
ferro and antiferro spin bonds as found in the classical spin 
model and density functional studies. 



FIG. 2: Left: nearest neighbor spin correlation function as 
calculated from the Monte Carlo simulation, the spin-orbital 
decoupling scheme, and the classical model. The correlation 
functions are normalized to their zero temperature values (- 
3/8, -1/4 and -1/8, respectively). In the first two cases spin 
valence bonds are formed, in the last case small residual inter- 
actions due to Hund's rule coupling will cause the formation 
of a chain of alternating ferro and antiferro bonds, see Fig. 3. 
Right: fit of the Monte Carlo results to the experimental spin 
susceptibility. 

so that the decoupled part of Hst — H^ c + H® c with 
fl|, = 47 £ (« + (-!)'<&) Si Si 

(i,3) 

+ ^(4J 8s + J CF )T t z . (5) 
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We focus first on the orbital part of this decoupled 
Hamiltonian. Since the only operator appearing in H® e is 
T z , it is effectively a classical Hamiltonian, closely related 
to the Ising model. This makes it possible to construct 
a transfer matrix that makes it possible to obtain an ex- 
act expression for the partition function. From there we 
can obtain any other thermodynamic quantity related to 
the orbital physics of H® e - The spin part of the decou- 
pled Hamiltonian (JSJ) cannot be diagonalized exactly as it 
stands. Therefore, following a canonical approach to the 
spin-Peierls Hamiltonian |l9j, we consider only the XY 
part of the spins. Then we can use the Jordan- Wigner 
transformations to diagonalize the remaining part of H^ e , 
and find all spin related thermodynamic quantities as 
modeled by H° c . 

Finally, we have considered the full Hamiltonian (J3J) 
without quantum spin fluctuations, i.e. taking only into 
account the z components of the spin operators. The mo- 
tivation for considering this limit is that in bandstructure 
calculations quantum fluctuations of this kind are effec- 



tively discarded. Since this Hamiltonian is classical, we 
can construct the transfer matrix that will allow us to 
find the partition function and thus all thermodynamic 
quantities of this system. 

Results. In Fig. 1 the spin and orbital ordering temper- 
atures are shown as a function of the crystal field splitting 
Jcf, as deduced from their respective susceptibilities in 
the Monte Carlo calculations. The transition tempera- 
tures for magnetic ordering coincide with the ones for or- 
bital ordering, as long as the crystal field splitting is not 
too large, i.e., when Jcf/ J ^ 1- The calculations show 
that the orbital and magnetic ordering occur at different 
temperatures only for a crystal field that is larger than 
the exchange coupling, i.e. when Jcf/ J ^ 1. In that 
case T° is proportional to Jcf, whereas T s tends to a 
constant. The same conclusion we find from decoupling 
calculations, shown in the right pannel of Fig. 1. 

We now invert the argument above. Experimentally 
the drop in spin susceptibility and the structural changes 
that are indicative of orbital ordering occur both at 
Too — 210 K. As experiment shows that these tem- 
peratures coincide, the splitting of ti g levels can thus not 
be larger than the superexchange coupling of neighbor- 
ing spins: evidently titanium pyroxene is in the regime 
where Jcf I J < 1. This conclusion can be further quan- 
tified. In Fig. 2, right pannel, we compare the intrinsic 
experimental spin susceptibility from Ref. |l(J and our 
Monte Carlo results. From fitting the experimental sus- 
ceptibility we find Jcf/ J = 0.7 ± 0.1, J = 156K and a 
gyromagnetic moment g — 1.85. 

We also calculated the nearest neighbor spin and or- 
bital correlation function C s jC° as a function of tem- 
perature with all three methods. The results for the spin 
correlator are shown in figure Fig. 2, left pannel. The 
corresponding orbital correlators C° indicate that in all 
cases there is a uniform and parallel ordering of the or- 
bitals at low temperatures and are not shown. The corre- 
lators in Fig. 2A are, for clarity, normalized to their zero 
temperature value. Clearly spins tend to align antiparal- 
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lei on neighboring bonds and form a valence bond solid, 
as shown in Fig. 3A. The Monte Carlo results show that 
this tendency is very strong and that even above the or- 
dering temperature spin dimer correlations are substan- 
tial. It is interesting to note that when the spin and 
orbital fluctuations are decoupled, the spin correlation 
function C s is appreciably reduced at finite tempera- 
tures. Thus the simultaneous fluctuations of spin and 
orbital degrees of freedom actually strengthen the spin 
valence bond solid. This seems counter intuitive, but has 
a straightforward physical cause. A spin flip in the va- 
lence bond state transforms a singlet dimer into a triplet, 
which is a strong ferromagnetic bond that causes an large 
reduction of the average antiferromagnetic spin correla- 
tions. Suppose, however, that the flipped spin can at the 
same time change orbital. After such a combined spin- 
orbital excitation, the spins that originally formed a firm 
singlet valence bond, are effectively not interacting any- 
more. They have a vanishing contribution to the near- 
est neighbor spin correlation function because of their 
respective orbital configurations and thus are less dam- 
aging to the magnetic ordering than fluctuations that 
exclusively involve the spin degree of freedom. 

The importance of fluctuations becomes even more 
striking when we compare the classical to the quantum 
results. When we only keep in the Hamiltonian H$t only 
the z projection of the spin operators, as is the case in the 
density functional calculations, the spin correlation func- 
tion C s is drastically reduced at any finite temperature, 
see Fig. 2A. Only at the lowest temperatures long range 
spin/orbital ordering sets in, which indicates that the sys- 
tem is strongly fluctuating due to the presence of many 
low energy states that are quasi degenerate. In this case 
we expect that the system is very susceptible to any resid- 
ual spin interaction. Indeed, we find that if we include a 
finite Hund's rule coupling, which leads to a very small 
ferromagnetic interaction between neighboring spins that 
are in different orbitals, the spins immediately order fer- 
romagnetically along the 7r-like bonds, whereas the spins 
on the a- like bonds, with orbitals pointing towards each 
other, order antifcrromagnetically, see Fig. 3B. In this 
case nearest neighbor correlation function C s , which is 
now the sum of antiferro and ferro correlations, goes to 
zero at T = 0. The resulting ground state spin and 
orbital configuration (Fig. 3B) exactly matches the one 
found by Popovic et al. in Ref. [ljj on the basis of 
bandstructure calculations. We find, however, that the 
ferro bonds are weak and the antiferro ones are strong, in 
accordance with recent LDA+U calculations 01 ■ There- 
fore even in this classical limit the system cannot be re- 
garded to be a spin one chain. 

Conclusions. Results from Monte Carlo simulations, 
from a spin-orbital decoupling and a classical approxi- 
mation scheme of a microscopic Hamiltonian which in- 
corporates both spins and orbitals of the titanium chains 



in NaTiSi20g place a firm upper bound on the size of 
the crystal field splitting. This proves that in this system 
the orbital degrees of freedom are active at low energies 
and temperatures. The interacting spin-orbital Hamilto- 
nian describes the experimental magnetic suscpetibility 
very well. The calculated correlation functions show that 
the orbital-Peierls state is realized below Too = 210-ftT 
and that quantum fluctuations in both the spin and or- 
bital sector are essential for the stabilization of this spin- 
orbital solid. 
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